Fluctuations of an Atomic Ledge Bordering a Crystalline Facet 
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When a high symmetry facet joins the rounded part of a crystal, the step line density vanishes 
as yfr with r denoting the distance from the facet edge. This means that the ledge bordering the 
facet has a lot of space to meander as caused by thermal activation. We investigate the statistical 
properties of the border ledge fluctuations. In the scaling regime they turn out to be non-Gaussian 
and related to the edge statistics of GUE multi-matrix models. 
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Equilibrium crystal shapes typically consist of various 
flat facets connected by rounded surfaces. For a micro- 
scopically flat facet there must be an atomic ledge border- 
ing the facet. This border step could be blurred because 
of thermal excitations, but is clearly visible at sufficiently 
low temperatures |l|, |2j, |3| ■ While in the interior of the 
rounded piece of the crystal the step line density is of or- 
der one on the scale of the lattice constant, it decays to 
zero as the edge of a high symmetry facet is approached. 
If r denotes the distance from the facet edge, according 
to Pokrovsky-Talapov |j] the step line density vanishes 
as y/r. Thus there is a lot of space for the border ledge to 
meander, in sharp contrast to steps in the rounded part 
which are so confined by their neighbors that they fluc- 
tuate only logarithmically [5j . The goal of our letter is to 
explore quantitatively the statistics of border ledge fluc- 
tuations. To illustrate our set-up we display in Fig. ^ a 
typical configuration from a statistical mechanics model 
which will be discussed below. One clearly recognizes 
the three facets as joined through a single rounded piece. 
Our interest is the statistics of the uppermost ledge. 

Experimentally ledge fluctuations are an elegant tool 
to determine step energies 0, : One carefully prepares 
an island, linear size L and single atom height, on a high 
symmetry facet. Alternatively, one sputters an undercut 
island. The ledge bordering the island is well described 
by a random walk, which implies fluctuations of size yL 
7\. In contrast, as can be clearly observed from Fig. ^ 
the border ledge of a facet interacts with its neighbors 
and a random walk model is not appropriate. In fact as 
our main result we will establish that the border ledge has 
fluctuations of size L 1 / 3 with a non- Gaussian statistics. 

To gain some understanding of the origin of such 
anomalous fluctuations let us consider the terrace-ledge- 
kink (TLK) model, which serves as an accurate descrip- 
tion of a vicinal surface, i.e. a crystal cut at a small angle 
relative to a high symmetry crystal plane. The surface is 
made up of an array of ledges which on the average run in 
parallel and are separated by terraces. The ledges are not 
perfectly straight and meander through kink excitations, 
only constrained not to touch a neighboring ledge. One 
can think of these ledges also as discrete random walks 




FIG. 1: Crystal corner viewed from the [111] -direction. 

constrained not to cross, i.e. with a purely entropic re- 
pulsion. Such a line ensemble is very closely related to 
Dyson's Brownian motion, in which the random walks 
are replaced by continuum Brownian motions. As dis- 
cussed in 0, Q , the location of the steps at fixed random 
walk time t have the same distribution as the eigenval- 
ues of a GUE ((3 = 2) random matrix. On this basis it 
is expected that the ledge-ledge distance is governed by 
the GUE level spacing [lCf. This prediction is verified 
experimentally [111] , however with deviations from (3 — 2 
which are attributed to long range elastic forces mediated 
through the bulk of the crystal and not included in the 
TLK model. 

If in the TLK model one retains the lattice structure 
in the transverse direction and makes the continuum ap- 
proximation in the direction along the ledges, then the 
ledges can be regarded as the world lines of free fermions 
in space-time ZxM UM- The world lines are piecewise 
constant and have jumps of only one lattice spacing. 
Consequently the transfer matrix has a nearest neigh- 
bor hopping term and the Pauli exclusion principle guar- 
antees entropic repulsion in the sense that ledges never 
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FIG. 2: Top lines for a TLK model with volume constraint. 
cross. 

The TLK model, in the version as just explained, has 
no facet. The crystalline surface has a constant average 
slope. Slope variations can be enforced through a volume 
constraint. For this purpose we introduce the "occupa- 
tion" variables rjj(t), \j\ < N, \t\ < T, in the surface 
patch [-N,-N+ 1,...,N] x [— T, T]: r/j(t) = 1 if there 
is some ledge passing through (j,t), and r]j(t) = oth- 
erwise. In these variables, up to an overall constant, the 
crystal volume is given by 



A v = dt JVj(t) 

J —T - ?vr 



(1) 



j=-N 



and volume constraint means to have an ensemble of 
ledges where the action A v is kept fixed. 

Without volume constraint the transfer matrix is gen- 
erated by a free fermion Hamiltonian with nearest neigh- 
bor hopping |l2^ . Imposing the volume constraint grand- 
canonically adds to the fermionic action the term \~ 1 A V 
with a suitable Lagrange multiplier A -1 . Thereby the 
nearest neighbor hopping Hamiltonian is modified to 



H F = 
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ajoj+i — a] +1 aj + 2a]ttj + ~aj<Zj). (2) 



a,j, resp. at, is the annihilation, resp. creation, operator 
at lattice site j 6 They satisfy the anticommutation 
relations {a.;,aj} = 6y , {<2i,a,} = = {aj,aj}. In (J2J 
we have taken already the limit N — *■ oo. The transfer 
matrix is e~ tHp , t > 0, and in the limit T — > oo one 
has to compute the ground state expectations for Hp. A 
macroscopic facet emerges as A — > oo. 

In Fig. |3 we display a typical ledge configuration for 
the TLK model with volume constraint. There is no fur- 
ther ledge above the one shown and for j : ^ — oo ledges 
are perfectly flat and densely packed. 

Since a ledge corresponds to a fermionic world line, the 
average step density (r)j(t))\ = p\(j) is independent of 
t and given by a = (<*jQj)\ with (-)\ on the right 

denoting the ground state expectation for Hp. By the 
linear potential in steps are suppressed for large j. 
Hence the average surface height hj(t) at (j,t), relative 
to the high symmetry plane, equals 



(3) 



(a^aj)x can be computed in terms of the Bessel func- 
tion Jj(z) of integer order j and its derivative Lj(z) — 
(d/dj)Jj(t) with the result 

Px{j) = {a) aj ) x = A(L i _ 1+a[A] (2A)J J+2[A] (2A) (4) 

-Lj +2 [x] (2A)Jj_ 1+2 [a] (2A)) 

where [•] denotes the integer part. For large A the height 
hj(t) is of order A. Therefore we rescale the lattice spac- 
ing by 1/A. Then lim^oo A _1 /i[A r ] (At) = h cq (r,t) with 
the macroscopic equilibrium crystal shape 

for r < —2 , 



h eq (r- 2,t) = { 



i (r arccos(r/2) — \/4 — r 2 ) 
for -2 < r < 2 , 



(5) 







for r > 2 . 



Thus under volume constraint the TLK model has two 
facets, one with slope 1, the other one with slope 0, joined 
by a rounded piece. The upper facet edge is located at 
r = 0. It has zero curvature. Expanding near r = 
results in h cq (r,t) = — ^jr(— r) 3 / 2 , consistent with the 
Pokrovsky-Talapov law. 

With the exact result (@J it becomes possible to refine 
the resolution. The appropriate step size is A 1 / 3 lattice 



constants. For the step density p\(j) = (a^a^x 
r = one finds 



lim A 1 / 3 p A (A 1 / 3 x) 

A — >co 



-xM{xf + Ai'(xf 



close to 



(6) 



Ai the Airy function. © has the asymptotics 
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for x 



-00 , 



(7) 



— exp(-4a; 3/2 /3) for 

&1TX 



Our real interest are the border ledge fluctuations. 
Clearly the border ledge is the top fermionic world line 
which we denote by b\(t). b\(t) takes integer values 
and is piecewise constant with unit size kinks. Since, 
at fixed t, the steps in the bulk have approximately the 
same statistics as a GUE random matrix, one would ex- 
pect that the transverse fluctuations of the border ledge 
equal those of the largest eigenvalue. Indeed, using the 
fermionic transfer matrix combined with an asymptotic 
analysis ^j|> one finds that 

lim Prob({6 A (0) < A 1/3 a}) = F 2 (a), a e R. (8) 

A — >oo 

In the random matrix community ^(d) is known as 
the Tracy- Widom distribution 14] . The correspond- 
ing probability density dF2(a)/da has an upper tail as 
exp(— |a 3 / 2 ) and a lower tail as exp(— jL|a| 3 ). 

In our context an experimentally more accessible quan- 
tity is the ledge wandering {[b\(t) — &a(0)] 2 ). In the limit 
of large A it has been computed in |l|J with the result 



k=j 



Var(6 A (t)-6A(0)) -A 2 / 3 .g(A- 2 / 3 i), 



(9) 



3 



using the short-hand Var(X) = ((X - (X)) 2 }. Thus the 
transverse fluctuations are on the scale A 2 / 3 . For small s 
the scaling function g(s) is linear in s, g(s) ~ 2|s|, indi- 
cating that for small, on the scale A 2 / 3 , separations the 
border ledge has random walk statistics. On the other 
hand g(s) saturates for large s, g(s) ~ 5(00) — c/s 2 , re- 
flecting that the border ledge fluctuations are station- 
ary (on the scale A 2 / 3 ). For the leading term one finds 
g(oo) = limA^oc A- 2 / 3 2(fe A (0) 2 ) = 1.6264. The sublead- 
ing coefficient c has recently been derived in 0, 0] with 
the result c — 2. 

Within the volume-constrained TLK model we arrived 
at an interesting prediction for the border ledge fluctua- 
tions. To be convincing we have to check against a more 
realistic model, for which we take the three-dimensional 
Ising model at low temperatures. At fixed crystal volume 
the equilibrium shape is then a cube with rounded cor- 
ners. Taken literally this model is still too complicated 
and we simplify through an SOS-type approximation by 
allowing only atomic configurations which are lattice con- 
vex. This means that, when the crystal is cut along any 
line parallel to the major axes, then the atoms fill a sin- 
gle interval (no holes). We use translation invariance 
to choose our coordinate system in such a way that the 
crystal lies in the positive octant of Z 3 with three of the 
facets coinciding with parts of the planes spanned by the 
three coordinate axes. If one now restricts attention to 
the piece of the crystal close to the origin, then the ac- 
tual crystal shape can be represented by a height function 
h(i,j), where i, resp. j, refers to the (100), resp. (010), 
axis. By construction 

(i) : h(i,j)>0, 

(ii) : h(i,j)>h(i + l,j), h(i,j)>h(i,j + l). 

The number of atoms missing, relative to the perfect 
cube, is 

V(h)= (10) 

i,j>0 

Thus the volume constraint translates into 
(Hi): V(h) — const. 

Every atomic configuration satisfying (i) — (iii) has the 
same number of broken bonds and thus the same en- 
ergy. Therefore our simplified version of the 3D Ising 
equilibrium droplet is to allow only atomic configurations 
which have a height function satisfying (i) — (iii) and to 
give them equal statistical weight. Note that our model 
is purely entropic. Fig. shows a typical sample with 
V(h) = 3x 10 5 . 

By projecting along the (111) direction the Ising cor- 
ner model is equivalent to tilings of the plane with rhombi 
of three distinct orientations. In this version the surface 
tension is computed in 17]. According to the Andreev 
construction the Legendre transform of the surface free 
energy yields the equilibrium crystal shape HE!. Con- 
venient formulas are available in |20| . where it is also 
established that, for the constraint V(h) = N, in the 



limit N — > 00 with corresponding lattice spacing A -1 / 3 
the equilibrium crystal shape is attained with probabil- 
ity one. From the implicit formula for the shape it can 
be deduced that near facet edges the Pokrovsky-Talapov 
law holds. The facet edge can be computed explicitly. If 
we consider the facet which lies in the (001) plane and 
denote the coordinate along (100) by r and the one along 
(010) by 6oo, then the macroscopic facet edge is given by 
600(7-) = -ln(l-e- T ), r>0. 

As observed in j^ij, the 3D Ising corner can be ana- 
lyzed through fcrmionic techniques. In particular, one 
can study the border ledge fluctuations. The details are 
rather intricate and given elsewhere [22^. Here we only 
report on those results which allow us to gain some under- 
standing of the universal properties of ledge fluctuations. 
We introduce the scaling parameter I by N = jC(3)^ 3 , 
C(3) = 1.202... being Apery's constant. The atomic 
border ledge position in the (001) plane is given by be(x), 
x = 0, 1, 2, ... . bi(x) takes positive integer values and 
is decreasing as bi(x + 1) < bi(x). The ledge has only 
South and East turns and meanders close to its asymp- 
totic mean £boo(x/£). We zoom at the fixed macroscopic 
edge point (It, ^600 (t)), t > 0. Upon proper rescaling 
|22|. one recovers exactly the same statistics as in JSJ. 
More precisely, for large £, 

Vax(bt(tT+x)-b t (tT)) = {\M) 2/3 g(-^ m x)> (11) 

with A = b'^r). (fTTfl differs from (JSJ in two re- 
spects. Firstly, to obtain the border ledge fluctuations 
one has to subtract the systematic mean. Since for our 
particular model the macroscopic facet edge is explicit, 
the subtraction is (be(£r + x)) — (be(£r)) = b'^^x + 
^b'^ (T)£~ 1 x 2 with negligible higher order corrections. 
Secondly, model-dependent properties enter indirectly 
through the coefficient A. Since g(s) = 2\s\ for small 
\s\, V&r(b e (£T +x)- bi(lr)) = A\ x I, A can be identified 
with the local wandering, resp. diffusion, coefficient. 

The border ledge of the TLK model and the 3D Ising 
corner have the same scaling behavior, which suggests the 
scaling to hold in greater generality. To obtain the form 
which properly distinguishes between model-dependent 
and universal properties we have to rely on a few notions 
from the thermodynamics of equilibrium crystal shapes 
[23| . Let us denote by h(x, y) the height of a vicinal sur- 
face relative to the high symmetry reference plane. We 
find it convenient to measure h in number of atomic lay- 
ers, whereas x, y are measured in a suitable macroscopic 
unit. Thus h is dimensionless and x,y have the dimen- 
sion [length]. Further let ^T/(u) be the surface free 
energy per unit projected area depending on the local 
slope u = Vh. Below the roughening transition / has a 
cone at u = and for small u behaves as 

/(u)- 7 (0)|u|+B(0)|u| 3 (12) 
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with 9 the polar angle of u |24J . The line stiffness 7 i s de- 
fined through j(9) — j(9) + j"(6). As argued in [25j, for 
short range surface models the Gaussian curvature of the 
equilibrium crystal shape has a universal jump across the 
facet edge, which implies the relation j(9)B(9) = tt 2 /6. 
Let us denote by / the Legendre transform of /. If 
J dxdy f '(V 'h(x , y)) is minimized under the constraint of 
fixed volume, then the resulting equilibrium surface is 
given by h(x,y) = £f(£~ 1 x,£~ 1 y), where I is the La- 
grange multiplier adjusted so to give the correct volume. 
h is convex downwards and has a convex facet lying in 
the x-y plane. The facet boundary is determined by "f(9) 
alone. Close to the facet edge, h = — |7ptg? 3 / 2 with 
d the normal distance to the facet edge, which defines 
the Pokrovsky-Talapov coefficient 7pt- Under Legendre 
transformation the angle 9 becomes the angle between 
the x-axis and the outer normal to the facet and, corre- 
spondingly, 7pt, the local curvature K, and the distance 
r of a point on the edge to the origin are parametrized 
through this angle 9. The relationship between 7 and B 
implies 

7 pt k = 2£- 2 n- 2 . (13) 

We return to the border ledge fluctuations close to a 
given angle 9q. For this purpose it is convenient to cen- 
ter the x-y axis coordinate system at t(9q) with the x- 
axis tangential and the y-axis along the inner normal to 
the facet. In this frame, we denote by y — b(x) the 
fluctuating border step. Then (b(x)) — ^k(9 )x 2 , in 
approximation. For sufficiently small \x\, still large on 
the scale of the lattice, b{x) is like a random walk and 
Var(6(ir) — 6(0)) = <7 2 |a;|, which defines the local wander- 
ing coefficient a 2 . Following [2^ it is natural to equate 
cr{9) 2 with the inverse stiffness j(9)~ 1 . This implies 

a 2 = k£, k = 7r 2 7 PT cr 4 /2 (14) 

valid for any point on the facet edge. 

The general scaling form is obtained now by using the 
TLK model as benchmark. Locally the border ledge per- 
forms a random walk with nearest neighbor hopping rate 
1, see ©, thus a 2 = 2. From © the PT coefficient is 
7pt = I/ttVI in our units. Using these two as model- 
dependent parameters yields the scaling form 

Var(6(x) - 6(0)) = (7r7 PT )- 4 / 3 5 ((7r7 PT ) 4 / 3 a 2 x/2). 

(15) 

Of course, through l|13|) . (|14l) . any other pair of model- 
dependent parameter can be used to reexpress l|15fl . 

As a control check, the Ising corner must also 
satisfy (|1 B|) . This is indeed the case with coeffi- 
cients k = £-V 2 , a 2 = &£,(t)(1 + 6 , oc (t) 2 )- 3 / 2 , and 

7PT = 2 1 /2^( T )-l^l/2 7r -l (1 + ^ (T )2)3/4. 

To summarize, the border ledge of a facet has fluctua- 
tions of size I 1 / 3 , thus much reduced in comparison with 



a simple random walk. We claim that the scaling form 
(|15fl is universal within the class of surface models with 
short range interactions. The scaling function g can be 
expressed through determinants of infinite dimensional 
matrices. Short and long distance behavior is known ex- 
plicitly. In H15|l there are two material parameters. Once 
they are determined experimentally, the functional form 
of the variance for the ledge fluctuations follows. 
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